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Abstract. Under general conditions, the equation g(x 1 , . . . ,x q ,y) = 
implicitly defines y locally as a function of x 1 , . . . , x q . In this article, 
we express divided differences of y in terms of divided differences of g, 
generalizing a recent formula for the case where y is univariate. The 
formula involves a sum over a combinatorial structure whose elements 
can be viewed either as polygonal partitions or as plane trees. Through 
this connection we prove as a corollary a formula for derivatives of y in 
terms of derivatives of g. 

1. Introduction 

Divided differences can be viewed as a discrete analogue of derivatives 
and are commonly used in approximation theory, see [1] for a survey. 

Recently, Floater and Lyche introduced a multivariate chain rule for di- 
vided differences [5], analogous to a multivariate form of Faa di Bruno's 
formula for derivatives [3,4,7]. In Theorem 1 in [8], this chain rule was 
applied to find an expression for divided differences of univariate implicit 
functions, thereby generalizing a formula by Floater and Lyche for divided 
differences of the inverse of a function [6] . 

In Theorem 3, the Main Theorem of this paper, we generalize Theorem 1 
in [8] to divided differences of multivariate implicit functions. More precisely, 
for some open box U C R q and open interval V C R, let y : U — > V be a 
function that is implicitly defined by a function g : U x V — > R via 

(1) 5 (x,y(x))=0, ||(x,y(x))/0 V x € U. 

Then the Main Theorem states that, for any rectangular grid 
{xq, . . . , x\ x } x • • • x {xq, . . . , x q nq ] C U, 

we can express the divided difference [xj, . . . , x 1 ^ ; • • • ; Xq, . . . , x«J y as a 
sum of terms involving the divided differences of g. 

In the next section, we define these divided differences and explain our no- 
tation. In Section 3, we apply the multivariate chain rule to derive a formula 
that recursively expresses divided differences of y in terms of divided differ- 
ences of g and lower-order divided differences of y. This recursive formula 
is "solved" in Section 4, in the sense that it is used to derive a closed- form 
nonrecursive formula that expresses divided differences of y solely in terms 
of divided differences of g. This formula is stated in the Main Theorem as a 
sum over polygonal partitions. It is shown in Section 5, that such polygonal 
partitions correspond to plane trees of a certain type, giving rise to an alter- 
native form of the Main Theorem. Switching between these combinatorial 
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structures, we are able to prove as a special case in Section 6 a generalization 
of a formula by Comtet, Fiolet, and Wilde for the derivatives of y in terms 
of the derivatives of g. 

2. Divided differences 
Consider a function y : U — > M. defined on some open box 

(2) U= (ai,6i) x ••• x (a q ,b q ) cR q . 

Suppose that, for some integers n±, . . . , n q > and all j = 1, . . . , q, we are 
given points Xq, . . . , x 3 nj 6 (aj, bj) satisfying aj < x < • • • < x 3 nj < bj. The 
Cartesian product 

(3) {xq, • • • , x\ x } x • • • x {xq, . . . , x q nq ] 

defines a [rectangular] grid of points in U. The divided difference of y at 
this grid, denoted by 

(4) [xj, . . . , x l ni ; • • • ; Xq, . . . , x^J y, 

can be defined recursively as follows. If n± = ■ ■ ■ = n q = 0, the grid consists 
of only one point (xj, . . . , Xq), and we define \x\; • • • ; xf^y := y(x\, . . . , Xq) 
as the value of y at this point. In case nj > for some 1 < j < q, we can 
define (4) recursively by 

(5) — ^o) [ X 0i • • • 5 x ni i " " " i x 0> • • • > x n 9 ] 2/ = 

[q™ 1 q^I • . . . • <T»*^ f y^ ry3 nr-3 • . « « • rr*i nr-Q 1 ri i 

Aq, . . . , x ni , , j,q, x l5 . . . , ^ ni _\i -i- nj , ) x 0' ■ • • ' ^nqj y 

r l l i i j j i a l 

'y> • • /ya*' /ya*' ryw /yW * • /yt J /y» ^ i 

|X , ■ ■ ■ , X m , , Xq, Xj, . . . , X n __ 1 , Xrij, , Xq, . . . , X„J (/, 

where the hat signifies omission of a symbol. If several of the nj are greater 
than zero, the divided difference (4) is uniquely defined by any of these 
recursive formulas. We refer to the dimensions (m, . . . ,n q ) of the grid as 
the order of the divided difference in Equation 4. 

For any a = (cti, . . . , a q ), b = (b±, . . . ,b q ) £ N q , write a < b whenever 
o-j < bj for every 1 < j < q. Additionally, we write a < b whenever a < b 
and a / b. In this manner, the symbol < defines a partial order on N 9 . We 
use the notation 

[x : a, h]y := \_x\ v x\ 1+1 , . . . , x^; ••• ; x^, x^ +1 , . . . , x q J y 

for the divided difference of y with respect to the grid of all points with 
indices "between a and b" . 

Divided differences of the function g : U x V — > M. in Equation 1 are 
defined similarly. For these functions, however, we stress the distinction 
between the variables x , . . . ,x q and the variable y by replacing the final 
semi-colon by a bar in our notation. 

As the notation of Equation 4 quickly grows cumbersome, we shall more 
often than not shorten the notation for divided differences to one that in- 
volves just the indices, 

(6) 

[k,o---k, sl ; ■■■ ;i q ,o---i q , Sq ]y-= [4, > ■ ■ • >4, S1 ; ■■■ ;< ' • ■ ■ > x LJf' 
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(7) [h,o ■ ■ ■ h, 31 ; • • • ; i q ,o ■ ■ ■ i q , Sq \jo ■ ■ ■ it] 9 ■= 

r i i . . q i I i 

L^ii.o' • • • ' X ii, si ' "' ' X i q ,o } ' " ' ' X i q ,s q \yjoi ■ ■ ■ ' Vjt\9- 

We can let some of the points coalesce by taking limits, as long as y is 
sufficiently smooth. In particular, letting all points in the grid coalesce to a 
single point xo = (xj, . . . , Xq) yields, for any tuple n = (m, . . . , n q ) € N 9 , 

r i i <? 9i _ 1 dH y i \ 

[x , . . . ,Xq] ■■■ ;x ,...,x \y— -j^-^rl x oJ- 

^^^^^^^ ^ ^ S 11 • C-A. 

ni+1 n q +l 

Here the derivatives are written in multi-index notation, |n| := rt\ + ■ ■ ■ + n q , 
and n! := n\\ ■ ■ ■ n q \ . Letting, in addition, the y- values coalesce to a single 
point yo yields 

[xq^Xo; ••• ;4 1 _^4\yo 1 _ : ^yo]g= — g xI1 ^ m ( x o^o)- 

ni+l n,+l m+1 

3. A RECURSIVE FORMULA 

Let y and y be related as in Equation 1. In this section, we derive a for- 
mula that expresses divided differences of y recursively as divided differences 
of g and lower-order divided differences of y. 

Consider a composition of functions M. q — > U x V — ^> M, where we write 
f : x i — > (/ 1 (x), . . . , / 9 (x), /'(x)). Here /' does not denote the derivative 
of some function / but is simply notation for the last component of f . Let 
be given a nonzero tuple n = (m, . . . , n q ) 6R' and a grid of points 

{x; : < i < n} = {sj, . . . x ••• x {x q , . . . , x q nq } C U, 

where Xq < ■ ■ ■ < xP nj for j = 1, . . . , q. Note that we allow for these coor- 
dinates to coincide. Let /.' := /'(x;) and f? := /-'(x;) for j = 1, . . . , q and 
< i < n. From [5, Theorem 2], we have the following multivariate chain 
rule in case f and g are sufficiently smooth, 

(8) [x:0,n]G,of) = £ £ £ 

fc=l O=io<---<ifc=n 0=jo<---<jr 9 +i=fc 
[4o ' 4'o+i ' ■ ' • ' 4/i ' " ' ' -^-i ' ^,-1+1 ' • ' • ' ^ J 4, ' 4,+i ' • • • ' 4,+l ] 9 

(1 jr \ I 3q+l \ 

II x:i ; ,.i/.r n , 
r=lj=>_i+l / \j=j q +l J 

where an empty product is considered to be 1 . (The formula in [5] includes a 
term for k = 0, but this term doesn't show up because we assumed n ^ 0.) 
For any integer q > 1, let us define a path [in Z g ] as a finite sequence 
(po, pi, • • • , p n ) of points in Z q for which po < pi < • • • < p n - We found 
it helpful to think of = jo < • • • < j q +i = k as indices along the path 
= io < • • • < ifc = n; see Figure 1 . 

Next, let f : x = (x 1 , . . . , x q ) \ — > (x, j/(x)) define the graph of a function 
y that is implicitly defined by g as in Equation 1. Let {ei, . . . , e q } denote 
the standard basis of R q , and let 1 < j < k and 1 < r < q be as in 



1 



GEORG MUNTINGH 




Figure 1 . Any sequence = io < ii < • • • < ifc = n repre- 
sents a path from to n. 

Equation 8. It follows directly from Equation 5 that the divided difference 
[x : ij_i,ij]/ r ' of the coordinate function f r : (x , . . . , x q ) i — > x r is equal 
to 1 whenever ij — = e r , and zero otherwise. The only choices of 

= jo < • • • < jq+i = k that yield a nonzero term in Equation 8 are 
therefore those satisfying 

(9) ij - ij-i = e r , for j = j r _i + 1, . . . , j r , r = l,...,q. 

Alternatively let (si, . . . , s q , t) := (ji - jo,..., j q - j q -!,j q+1 - j q ) be the 
sequence of jumps in the sequence (jo, . . . , jq+i). In terms of these jumps, 
Equation 9 is equivalent to the statement that the path ig < • • • < ijt starts 
with 

si steps of ei, followed by 
S2 steps of e2, followed by 

s q steps of e q , followed by 
t arbitrary steps. 

Let us call any tuple (s±, . . . , s q , t) with this property compatible with (io, . . . , ifc), 
or simply compatible if it is clear which sequence (io, • • • , ife) is referred to. 
Note that such a tuple forms an integer partition k = s% + • • • + s q + t. 
Equation 8 thus implies 

M 

(10) [x:0,n]g(;y(-))=Y / E E 

k=l 0=io<---<ifc=n compatible 
(s 1 ,...,s q ,t) 

\s\+t 

[01---si; ••• ;01---Sg|i[ s [i[ s | +1 ---i| s | +t ]p [x.:ij-i,ij]y, 

j=\e\+l 

where we used the shorthand notation 

[01 ■••si; ••• ;01---s g |i| s |i| 8 | +1 ---i| s | +t ]5 

= [x , x 1: . . . , x Sl ; ■■■ ; Xq, x^, . . . , |yi| s | , yi |s|+1 , . . . , yi| s|+t ] 

Hi := j/(xi), from Equation 7, for the divided differences of g. 

If y is implicitly defined by g as in Equation 1, the left hand side of 
Equation 10 is zero. Suppose n has length |n| = 1. Then n = e r for some 

1 < r < q. In this case, the right hand side of Equation 10 comprises two 
terms with k = 1 and = io < ii = e r . One finds 

0=[x:0,e r ]g(;y(-)) = 
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[0; ••• ;0|0e r ] 5 [x:0,e r ]y+[0; ••• ;0;01;0; ••• ;0\e r ]g, 

r-1 

or, equivalently, 

r-1 

r n , ^~0;Qi;Q; ••• ;Q|er]g , 1 

(Rl) [x:0,e r]y = [Q; ...^Oe^ ' **r=l,...,q. 

For example, when y is a function of q = 2 variables, this equation represents 
the two formulas 

• ' 1 ' J " [0;0|0ei] 9 ' 1 " [0;0|0e 2 ] 9 

Now suppose n has length |n| > 1. There is only one term in the right 
hand side of Equation 10 with k = 1. This term is given by 

= io < ii = n, si = ■ ■ ■ = s q = 0, t = 1 

and involves the highest order divided difference of y. Isolating this divided 
difference yields a formula that recursively expresses divided differences of 
y in terms of divided differences of g and lower-order divided differences of 



(R2) [x:0,n] y = £ £ £ 

k=2 0=io<---<ifc=n compatible 
(si,...,s q ,t) 



( [01---«i; ■■■ ; 1 ■ • - a 9 1 - ■ ■ i|.| + i]g \ 'S' , , 
V [0;- ;0|On ]s j nj* ■ ¥-1. 

Let us simplify this formula. The product in Equation R2 can be split 
into two products 

|s|+t |s|+t 

(12) Yl ■ t x : V-i'^ly- 

j=|s|+l j=|s|+l 
lij-ij-iN 1 |ij-ij-i|>2 

From Equation Rl it follows that each divided difference in the first product 
can be expressed as a quotient of divided differences of g. Our ultimate goal 
is to express the left hand side of Equation R2 solely in terms of divided 
differences of g. To achieve this, it seems natural to split the right hand side 
into a part that can directly be expressed in terms of divided differences of 
g and a remaining part involving higher-order divided differences of y. The 
former part can be expressed by introducing, for every sequence 

(io, • • ■ , ii) = h < ■ ■ ■ , if) = h < ■ ■ ■ < (4> • • • , i\) = 
a symbol {io • • • ik}g for the expression 

[jo — (jo+_ai)j ••• ;ig---(ig + s g )|i| g |---i| 8 | +t ]g \ 



(13) E 



compatible 

(si,...,S q ,t) 
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lj — \ = Gr 

involving only divided differences of g. Whenever it is hard to separate visu- 
ally the multi-indices io, . ■ • , ife, we write {io, . • • , i/c}<? instead of {io • • • ifc}<7- 
The divided differences [x : ij_i,ij]y that appear in the second product 
of Equation 12 satisfy \ij — > 2. As Equation 9 guarantees that this 
cannot happen for j < j q = |s|, we might as well start the product of these 
remaining divided differences at j = 1 instead of at j = |s| + 1. This has the 
advantage of making the expression independent of |s|. Equation R2 can 
thus be written in the concise form 

|n| ^ k 

(R2') [x:0,n]y = ^ ^ {i • • • i k }g ]J [x : i^, 

k=2 O=i <--<i fc =n j=l 

\ij-ij-l\>2 

Equation Rl gives a formula for [x : 0, n]y when |n| = 1. Let us consider 
Equation R2' for the case that |n| = 2. For such n, either n = 2e r with 

1 < r < q, or n = e r + e s with 1 < r < s < q. In the examples below we 
compute [x : 0, n]y for these two cases, assuming q = 2 to simplify notation. 

Example 1. Suppose n = 2ei (the case n = 2e2 is similar). The only 
possible path = io < • • • < i& = n with k = 2 in Equation R2' is given 
by < ei < 2ei. As for such a path the product in Equation R2' is empty, 
one has [x : 0, 2e\]y = {0, ei, 2ei}g. 

To compute {0, ei, 2ei}g, we need to find out which integer partitions 

2 = s\ + S2 + t are compatible with this path. These are precisely the triples 
(si, S2,t) for which 

(14) ei = ii - i = i 2 - ii = ■ ■ ■ = i si - i si -i, 

(15) e 2 = i si +i — i si = isi+2 — isi+i = • • • = i Sl + S2 — i ai + S2 _i, 

where the first (respectively second) statement is considered to be trivially 
satisfied whenever s\ = (respectively S2 = 0). As both ii — io and i2 — ii are 
equal to ei, the first condition is automatically satisfied. As neither ii— io nor 
'i2 — ii is equal to e2, necessarily S2 = 0. It follows that there are three triples 
(si,s 2 ,t) = (0,0,2), (1,0,1), (2,0,0) compatible with (0,ei,2ei). Each of 
these sequences corresponds to a term in {0, ei, 2e\}g, and we conclude that 

(16) [x:0,2ei]y = {0,ei,2ei}<7 = 

[0;0|0,ei,2ei]g [01;0|ei]g [12;0|2e 1 ]g 

[0;0|0,2ei] 5 [0; 0|0, e^g [1; 0|ei, 2e l )g 
[01;0|ei,2ei]g [12;0|2ei]g 

[0;0|0,2ei]g [1; 0|e 1; 2ei]g 
[012;0|2ei]g 
[0;0|0,2ei]g" 

Example 2. Suppose n = ei + e 2 . Equation R2' is a sum over the two 
possible paths < ei < ei + e 2 and < e 2 < e\ + e 2 . 

A triple (si, s 2 , t) is compatible with the path < ei < ei + e 2 precisely 
when Equations 14 and 15 hold. For this path, the first equation is equivalent 
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to si either being or 1. If si = 0, then the second equation implies that 
s 2 = 0. If si = 1, on the other hand, the second equation implies that S2 is 
either or 1. One finds three triples (si,s 2 ,t) = (0,0,2), (1,0,1), (1,1,0) 
compatible with (0, ei,ei +e2), yielding 

{0,ei,ei +e 2 }g = 

[0; 0|0, ei, ei + e 2 ]g [0 1; Ojgijff [1; l|ei + e 2 ]g 

[0;0|0,ei +e 2 ]g [0; 0|0, e^g [1; 0|ei, d + e 2 ]g 
[01;ei|ei,ei + e 2 ]g [1; l|ei + e 2 ]g 

[0;0|0,ei +e 2 ]g [l;0\e 1 ,e 1 + e 2 ]g 
[01;01|ei+e 2 ] g 
[0;0|0,ei + e 2 b' 

Similarly, a triple (si, s 2 ,t) is compatible with the path < e 2 < e\ + e 2 
precisely when Equations 14, 15 hold. For this path, however, the fact that 
i 2 — ii = ei comes after ii — io = e 2 implies that si = 0. One finds two 
triples (s\,s 2 ,t) = (0,0,2), (0,1,1) compatible with (0, e 2 , e\ + e 2 ), yielding 

{0,e 2 ,ei +e 2 }g = 

[0; 0|0, e 2 , ei + e 2 ]g [0; l|e 2 ]g [0 1; l|ei + e 2 ]g 
[0;0|0,ei +e 2 ]g [0; 0|0, e 2 ]g [0; l|e 2 , ei + e 2 ]g 
[0;01|e 2 ,ei + e 2 ]g [0 1; l| ei + e 2 ]g 
[0;0|0,ei +e 2 ]g [0; l|e 2 , ei + e 2 ]g' 
As for both paths the product in Equation R2' is empty, it follows that 

(17) [x : 0, ei + e 2 ]y = {0, ei, ei + e 2 }g + {0, e 2 , e x + e 2 }g. 

4. A FORMULA FOR DIVIDED DIFFERENCES OF IMPLICIT FUNCTIONS 

Let y be implicitly defined by g as in Equation 1. In this section we derive 
a formula that expresses divided differences of y solely in terms of divided 
differences of g. For n with |n| = 2, Equation R2' immediately yields the 
two formulas 

(18) [x:0,2e r ]y = {0,e r ,2e r }( ? , 

(19) [x : 0, e r + e s ]y = {0, e r , e r + e s }g + {0, e s , e r + e s }g, 

where 1 < r < s < q and the expressions {io • • • i&} are defined in Equation 
13. For n with |n| = 3, one can distinguish three cases: n = 3e r , n = 
2e r + e s , and n = e r + e s + e^, with 1 < r, s, t < q distinct. Let us compute 
[x : 0, n] for these n to get a feel for what a general formula should be. 
Repeatedly applying Equation R2' yields 

(20) [x : 0, 3e r ]y = 

{0, e r , 2e r , 3e r }g + {0, 2e r , 3e r }<7 • {0, e r , 2e r }g 
+{0, e r , 3e r }<7 • {e r , 2e r , 3e r }g 

(21) [x : 0, 2e r + e s ]y = 

{0, e r , 2e r , 2e r + e s }^ + {0, 2e r , 2e r + e s }g ■ {0, e r , 2e r }g 
+{0, e r , 2e r + e s }g • {e r , 2e r , 2e r + e s }g 

+{0, e r , e r + e s , 2e r + e s }g + {0, e r + e s , 2e r + e s }g ■ {0, e r , e r + e s }g 



8 



GEORG MUNTINGH 



+{0, e r , 2e r + e s }g • {e r , e r + e s , 2e r + e s }g 

+{0, e s , e r + e s , 2e r + e s }g + {0, e r + e s , 2e r + e s }# • {0, e s , e r + e s }g 
+{0, e s , 2e r + e s }g • {e s , e r + e s , 2e r + e s }# 

(22) [x:0,e r + e s + ei]y = 

{0, e r , e r +e s , e r +e s +eJ#+{0, e r +e s , e r +e s +e t }g-{0, e r , e r +e s }g 

+{0, e r , e r + e s + e t }# • {e r , e r + e s , e r + e s + e t }# 

+{0, e r , e r +e t , e r +e s +e t }5i+{0, e r +e t , e r +e s +e t }g-{0, e r , e r +e t }g 

+{0, e r , e r + e s + e t }g ■ {e r , e r + e t , e r + e s + e t }g 

+{0, e s , e r +e s , e r +e s +e t }g+{0, e r +e s ,e r +e s +e t }g-{0, e s , e r +e s }g 

+{0, e s , e r + e s + e t }g ■ {e s , e r + e s , e r + e s + e t }g 

+{0, e s , e s +e t , e r +e s +e t }g+{0, e s +e t , e r +e s +e t }5-{0, e s , e s +e t }g 

+{0, e s , e r + e s + e t }g ■ {e s , e s + e t , e r + e s + e t }g 

+{0, e t , e r +e t , e r +e s +e t }g+{0, e r +e t , e r +e s +e t }g-{0, e t , e r +e t }g 

+{0, e t , e r + e s + ejy • {e t , e r + e t , e r + e s + etjy 

+{0, et, e s +e t , e r +e s +e t }g+{0, e s +e t , e r +e s +e t }g-{0, e t , e s +e t }g 

+{0, e t , e r + e s + ejy • {e t , e s + e t , e r + e s + e t }g 

These three formulas exhibit a remarkable pattern. For every choice of the 
path = po < pi < P2 < P3 = n, we seem to be getting a sum 

{P0PlP2P3}5 + {PoP2P3}5' • {PoPiP2}5 ( + {poPiPalfi' • {P1P2P3I5- 

This expression bears a striking resemblance to the right hand side of the 
univariate formula 

[012% = {0123} 5 + {023}g {012}g + {0l3}g{l23}g 

established in [8, Theorem 4]. This suggests that, for general n = (m, . . . , n q ), 
the divided difference [x : 0,n]y is a sum of ( n ^~ " + ™ q ) univariate formulas, 
one for each choice of the path = po < • • • < P| n | = n. See Figure 2a for 
an example of such a path. 

Theorem 3 casts this suspicion into a precise form. In order to state this 
Theorem, we introduce some notation for polygon partitions. With a se- 
quence of labels po, Pi, • • ■ , p n we associate the ordered vertices of a convex 
polygon. A partition of a convex polygon is the result of connecting cer- 
tain pairs of nonadjacent vertices with straight line segments, none of which 
intersect. We denote the set of all partitions of the polygon with vertices 
p ,pi, • • • ,p n by P(po,Pi, • • • ,Pn)- The points p ,Pi, • • • , p n and line seg- 
ments (between either adjacent or nonadjacent vertices) form the vertices 
and edges of a plane graph. As such, every partition ir G V(po, pi, . . . , p n ) 
is described by its set F(tt) of [oriented] faces, which does not include the 
unbounded face. Each face / £ ^(vr) is represented by a subsequence 
/ = (vo, vi, . . . , Vfc) of the sequence (po, Pi, . . . , p n ) of length at least three. 
We let E(ir) denote the set of edges in ir, each of which is represented by 
a subsequence (vo, Vi) of (po, Pi, • • • , Pn 

) of length two. An edge (v ,Vi) 
is called an outer edge if the vertices Vo, Vi are adjacent. If Vo,Vi are not 
adjacent, (vo,Vi) is called an inner edge. Figure 2b depicts an example of 
such a partition of a convex polygon. 
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P 5 Pe P 7 =( n i. n 2) 
"p 4 

|Pi *P 2 P 3 

(0,0) = Po ' 

(A) (B) 

Figure 2. For n = (ni,n 2 ) = (4,3), the figure to 
the left depicts a choice of a path = po < • • • < 
P| n | = n. The figure to the right shows a parti- 
tion of the convex polygon corresponding to this path 
with faces (p , Pi, P2, P7), (P2, P3, P5, Pe, P7), (P3, P4, Ps), in- 
ner edges (p2, P7), (p3, ps) (drawn solid), and outer edges 
(Po,Pi), (P1,P2), • • • , (P6,P7), (P0,P7) (drawn dotted). 

Armed with this notation for partitions of convex polygons, we are now 
able to state the Main Theorem of this paper. 

Theorem 3 (Main Theorem). For n with |n| > 2, 
(23) 

[x:0,%= E E II { v ° Vl ' 

0=Po<Pi<"-<P|n|=n 7reP(po,pi,-,P|„|) (v ,vi,...,v r )eF(7r) 

To prove Theorem 3, our plan is to use Equation R2' recursively to express 
[x : 0, n]y solely in terms of divided differences of g. Before we proceed with 
this proof we assign some visual meaning to Equation R2' to highlight the 
backbone of this proof. We call a sequence i = (io, ii, . • • , ifc) a subpath of 
P = (Po, Pi, • • • , Pn) and p a superpath of i, whenever 

io = Pi < Pio+i <■■ <h = Ph < Ph+i < ■ ■ ■ < ifc = Pi k , 

for some increasing indices = Iq < l\ < • • • < = n. Every subpath i of p 
induces a partition in "P(po, Pi, • • • , Pn) whose set of faces comprises an inner 
face (io, ii, • • • , ifc) and outer faces ...,ij) = (pi^Pi^+x, ■■■,Pi J ) for 

every j = 1, . . . , k with |ij — > 2. See Figure 3b for an example. 

In general, a sequence (io, ii, • . . , ifc) has several superpaths (po, pi, . . . , p n ). 
Let us introduce some notation to consider simultaneously the partitions of 
the outer faces (each of which is a convex polygon itself) of i = (io, ii, . . . , ifc) 
for all these superpaths. We define 
fc 

V\ ■■= Yl II P(qo,qi,---,qm), 

i=l ij-_i=qo<— <qm=ij 

m:=\ij — ij_i|>2 

which represents a set of tuples of partitions, each entry in such a tuple corre- 
sponding to a partition of a path with steps in {ei, . . . , e q } from to ij for 
some j. For example, fori = (i ,ii, . . . ,i 5 ) = ((0,0), (0, 1), (1, 1), (2,2), (2,3), 
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(0,0) =i " 

(A) 



1 5 =(nj ,n 2 ) 



P3 P4 =i 3 




P0 =i P7 =i 5 

(B) 



Figure 3. For n = (ni,n 2 ) = (4,3), the figure to the 
left shows the points in the sequence i = (io, ii, . - - , is) = 
((0,0), (0,1), (1,1), (2, 2), (2, 3), (4, 3)), together with the 
paths traced out by its two superpaths p = (po, Pi, . . . , P7). 
The figure to the right shows the convex polygon correspond- 
ing to each of these paths. The sequence i gives rise to two 
outer faces (P2>P3>P4) and (ps,P6>P7)> which are drawn 
shaded in the figure. Depending on the choice of the su- 
perpath p, the former outer face is either (12, (1, 2), i 3 ) or 
(i 2 , (2, 1), is) , while the latter is equal to (i4, (3, 3), is) for 
both paths p. 



(4,3)), one has |ij — ij_i| > 2 only for j = 3,5 (see Figure 3a). There are 
two paths with steps in {ei, 62} from i2 = (1, 1) to i3 = (2, 2) and only one 
from 14 = (2, 3) to is = (4, 3). It follows that 

Vi = (V(i 2 , (1, 2), i 3 ) Up(i 2 , (2, 1), i 3 )) x V(h, (3, 3), i 5 ) 

= { ((i 2 , (1, 2), i 3 ), (i 4 , (3, 3), i 5 )) , ((i 2 , (2, 1), i 3 ), (i 4 , (3, 3), i 5 )) } . 

We now associate divided differences to these geometric objects. To each 
outer face (ij-i, ■ ■ ■ , ij) we associate the divided difference [x : ij-i, ij]y, and 
to each inner face (io, ii, • • • , ifc) we associate the expression {io • ■ ■ ifc}<?- For 
any sequence i that appears in the sum of Equation R2', the corresponding 
inner face therefore represents that part of Equation R2' that can be written 
solely in terms of divided differences of g, while the outer faces represent the 
part that is still expressed as a divided difference of y. 



Proof of Theorem 3. The proof is by induction on |n|. Equations 18-22 show 
that the formula holds for |n| = 2, 3. For a fixed |n| > 4, suppose the formula 
holds for all smaller |n| (but with |n| > 2). Consider the recursion formula 
from Equation R2'. As in each term k > 2, one has |L — ij— 1| < |n| for j = 
1, . . . , k. By induction, therefore, we can replace each divided difference [x : 
ij-i, ij]y in Equation R2' by an expression involving only divided differences 
of g. The symbol "Pi enables us to consider these expressions for [x : ij—i, ij}y 
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simultaneously, yielding 

k 

(24) H [x-.ij-i^y 

i=i 

|ij-ij-i|>2 

- LT £ £ n to-")* 

i =1 ij_i=qo<-<qm=ij 7r6'P(qo,...,q m ) (v ,...,v r )GF(7r) 

m—lij— ij_i|>2 

e n {vo---v r }5. 

(7ri,7r 2 ,...)GPi j>l 

(v„,...,v r )eF(7r,) 

For a given sequence = io < • • • < i& = n with k > 2, the set Pi can be 
identified with the set 

{vr G -P(p , • • • , P| n |) : = po < • • • < P| n | = n, i G F(tt)} 

by the bijection that maps any tuple (7Ti,7T2, . . .) in Pi to the partition ir 
with E(7r) = {i} U F(tti) U F(it2) U • • • . Applying this bijection to Equation 
24 and substituting the result into the recursive formula from Equation R2' 
yields 

[x : 0, n]y 

|n| 

= E E {io • • • ifc}5 E II { v o--- v r}5 

fc=2 0=i <-"<ife=n (7ri,7r 2 ,...)GPi j>l 

(vo,...,v r )eF(7r J ) 

= E E E n iv ...v r}9 

k=2 0=i <-<i fc =n 7r6-p(p ,...,P|„|) (v ,...,v r )eF(7r) 
a subpath of iGF(Tr) 
0=p <---<P| n |=n 

E E II {^■■-r}g. □ 

0=p <-<P|„|=n 7reP(po,..., P | n |) (v ,...,v r )eF(7r) 

5. Polygon partitions and plane trees 

While the compact nature of Equation 23 is useful to state and prove 
Theorem 3, it is less appropriate for finding a specialized formula for deriva- 
tives. In this section we adapt Equation 23 to a form better suited for this 
purpose. 

Following [9, p. 294], one defines a plane tree T recursively as a finite 
set of vertices such that one specially designated vertex is called the root, 
and the remaining vertices are put into an ordered partition (T±, . . . , T rn ) of 
m > disjoint nonempty sets Ti, . . . , T m , each of which is a plane tree. We 
recall the following Lemma, which appears as Proposition 6.2.1 in [10]. 

Lemma 4. For all integers m, n with m > n > 2, there is a bijection between 
the following two structures: 

• Plane trees with m vertices of which n are leaves, and all other ver- 
tices have at least two descendants. 
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• Partitions with m — n faces of a convex polygon with n + 1 vertices. 

We now explicitly describe this bijection. Suppose we are given a polygon 
partition in 7 ? (po,Pi, • • • , Pn) with m — n faces. To the edge (po,Pn) we 
associate a vertex that represents the root of our tree. As (po, p n ) is an outer 
edge, it belongs to a unique face. The other edges of this face are taken to be 
the descendants of the root vertex. As we are constructing a plane tree we 
need to order these descendants; the vertices correspond from left to right 
to the edges encountered when traversing the border of the face clockwise, 
starting at (po,Pn)- We then repeat this process for each of the new edges 
until we are out of edges. This construction yields a plane tree with n leaves 
corresponding to the outer edges (po, Pi), (pi, P2), ■ ■ ■ , (Pn-l> Pn) and m — n 
nonleaf vertices corresponding to the faces of the polygon partition. 

For example, for the partition in "P(po, Pi, • • • , P7) with set of faces 

{(P0,P5,P6,P7), (P0,P2,P5), (P0,Pl,P2), (P2,P4,P5), (P2,P3,P4)}, 

this bijection can be visualized as follows. 



P3 P4 P3 P4 Q O 




Here leaves are drawn as circles and nonleaf vertices are drawn as discs. 

Conversely, suppose we are given a plane tree with m vertices of which n 
are leaves, and all the other vertices have at least two descendants. Assign 
labels (po> Pl)j (Pi) P2)) • • ■ , (Pn-i, Pn) to the leaves as they are encoun- 
tered while traversing the tree depth- first from left to right. Recursively, we 
assign the label (pj , p» 1 , . . . , pi r ) to any nonleaf vertex whose descendants 
have labels of the form (p io , . . . ,p h ), (p h , . . . ,p i2 ), . . . , (pi^ , . . . , Pi r ) . The 
labels of the nonleaf vertices then coincide with the faces of a partition in 
V(po, pi, . . . , p n ), and the labels of the leaves correspond to the outer edges 
unequal to (po, p n ) of the full polygon. The following picture illustrates this 
construction with an example. 



(Po.Pl) (P1.P2) (Ps.Pe) (Pn.Py) 
OOOO OOOO 




Let 7~(po, Pi, • • ■ , Pn) denote the set of plane trees with n leaves for which 
all nonleaf vertices have at least two descendants. Represent each tree 
r G 7~(po, Pi, • • • , Pn) by its set V(r) of nonleaf vertices that are labeled as 
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above, and label the leaves correspondingly. Note that V(t) = F(tt) when- 
ever a tree r and polygon partition ir are related via the above bijection. 
From this it follows that Equation 23 can equivalently be stated in terms of 
plane trees as 

(23') 

[x:0,n]y= ^ ^ ]J {ioh ■ ■ ■ k} 9 ■ 

0=po<Pi<---<P| n |=n reT(po,pi,-.P|n|) (io,ii,-,ifc)eV(r) 

We wish to bring this equation into a form where we can distinguish 
the individual terms in the divided differences of g. For this, we replace 
T(po, • • • , P| n |) by a structure T'(po, • • • , P|n|) that encompasses all combi- 
nations of all terms in the expressions {io ■ ■ ■ ifc}<7- More precisely, every tree 
r in T(po, • • • , P| n |) wm be replaced by several trees, one for each term in 
Equation 13, each of which is an extension of r with an additional nonleaf 
vertex for every factor in the second line of Equation 13. 

Let r be one of the trees in T(po, • • • >P|n|) corresponding to a polygon 
partition it. Any nonleaf vertex v = (io, - - - , ifc) in V(r) defines, together 
with its direct descendants, a subtree * v of r called a [plane] star with 
root v, in which we remember which descendants were (non)leaves. Note 
that the bijection of Lemma 4 induces a bijection between F(n) and the 
set Stars(r) of stars of nonleaf vertices of r. A star * is said to be of type 
(s*, . . . , s*, t*) = (s*,t*), if the sequence of descendants of its root starts 
with 

s* leaves with labels (a, b) satisfying b — a = ei , followed by 
«2 leaves with labels (a, b) satisfying b — a = e2 , followed by 

s* leaves with labels (a, b) satisfying b — a = e q , followed by 
t* nonleaves. 

See Figure 4 for an example. Note that such a type does not exist for every 
star, as leaves can appear after nonleaves. 

For every integer partition k = s± + ■ ■ ■ + s q + t compatible with v = 
(io, . . . , ifc), we can extend * v to a tree t v '* by inserting an edge at the leaves 
among the final t descendants (i| s |, i| s |+i), (i|s|+i, i|s|+2), ■ ■ ■ , (i| s |+t-i, i|s|+t) of 
v. That is, we insert an edge for every factor in the second line of Equation 
13. Note that if there are no such factors, then tv'* = * v - Every star ★ in 
Tv* is then of some (necessarily unique) type (s*,i*). Using these notions, 
one can write 

{i ---ifc}5= Yj n 

compatible * e StarsfTv'*) 
(B,t) 

/ [#•••(# + *?); ••• ;^---(^ + <)|i| S r--W]^ 

where each star ★ is of type (s*, t*) = (s*, . . . , s*, t*) and has root (ig, . . . , i£*) , 
with 

(zj*, . . . ,Zq ) = ig < ■ ■ ■ ,i\ ) = ii < ■ ■ ■ < (ifc*> ■ ■ > = ifc*- 
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Figure 4. The figure to the left shows the path 
= (2e 1 ,3ei,4ei,4e 1 + e 2 ,5ei + 3e 2 ). The 
figure to the right shows a star * of type (sf , st,, t*) = (2, 1, 1) 
with root (ig, i*, i|, ij, i|), with the circles representing leaves 
and the discs nonleaves. 



Let T'(po) • • • jPlnl) be the set of plane trees obtained by taking a plane 
tree in T(po, • • • >P|n|) an d replacing each of its stars * v by r*' for some 
(s,i) compatible with v. Equivalently, T"'(po, • • • , Plnl) is the set of plane 
trees with leaves (po,Pi), (Pi)P2), (P|n|-i, Plnl) an d nonleaves labeled 
accordingly, for which each star * is of some type (s*, t*) 7^ (0, 1). Equation 
23' can then be stated as 

(23") [x:0,n]y= £ £ J] 

0=p <pi<---<P| n |=n r'eT'(po,Pi,.--,P|n|) *6Stars(r') 



1* . . . r„-i* _L o*V ... .. . (it* _i_ D *M5* 



' ' ' (*0 + s l^ ' ' ' > *0 



(*0 + fl $)l i fs*| •■• i fs*|+t*]5 



where again each star * is of type (s*, t*) = (sf, . . . , s*, t*) and has root 
•••>***)> with 

(<S* . . . ,<§*) = is < (<}*, • • ■ ,<f) = ^ < • • • < (<& • ■ • ,<$) = ij*. 

Example 5. To the path (po, Pi, P2, P3) = (0, ej.,ei + e2,2ei + e 2 ) corre- 
spond three trees in T(po, Pi, P2, P3), 



n 

(Pl ,P 2 ) (P 2 >P3 ) 



T-2 

(Po 'Pl ) (Pl >P2 ) 



73 



(P0,Pl)V/ . , V/ (P2,P 3 ) (Po.Pl) (Pl,P 2 ) (Pz.Ps) 

q HP1.P2.P3) (p >Pi>p 2 )V 



(Po Pl »P3 ) 




(P0.P2.P3) 




(P05Pl.P2.P3 



Let us consider the first tree t\. There are two tuples (si, t\) = (0, . . . , 0, 2), 
(s 2 , t 2 ) = (1) 0, . . . , 0, 1) compatible with the nonleaf vertex v = (po, Pi, P3), 
and we can extend * v to two corresponding trees 
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r, 

(Po >Pl 



Sl,tl 



(Po,Pi> 



(Pi ,P 2 >P3 I 




(P >Pl =P3 



_S2,*2 



(POfPl) (Pl>P 2 .P,3 




(P0.P1.P3 



Similarly, there are two tuples (s' l5 i^) = (0, . . . , 0, 2), (s' 2 , t' 2 ) = (0, 1, 0, ... , 0, 1) 
compatible with the other nonleaf vertex v' = (p1.P2.P3) of t±, and we find 
two trees 



s' v t[ 

V 

(Pl,P 2 ) (P 2 .P 3 ) 

9 9 



(P1.P2' 




(P 2 .P 3 ) 



(Pi ,P 2 ,P,3 



_ s 2'*2 
V 



(P 2 ,P 3 ) 

O 



(pi>p 2 ; 




(P 2 >P,3 ) 



(P1.P2.P3) 



corresponding to * v '. It follows that the tree T\ yields 2x2 = 4 different 
trees in T'(po, Pi, P2, P3)- Analogously, one can check that T2 yields 1x3 
trees and T3 yields 3 trees. 

We end this section with a Lemma that appears as Theorem 5.3.10 in [10] 
and is needed to compute the coefficients in Equation 25. 

Lemma 6. For any k £ N, there are 

1 fro + n H h r k 

ro + n-\ \-r k \ r 0) ri, . . . ,r k 

different plane trees with ro vertices with descendants (leaves), r\ vertices 
with 1 descendant, . . ., r k vertices with k descendants, and no vertices with 
more than k descendants. 



6. IMPLICIT HIGHER PARTIAL DERIVATIVES 

Whenever g and y are sufficiently smooth, coalescing the grid points in 
Equation 23" results in a formula for the derivatives of y in terms of the 
derivatives of g. We show that this formula generalizes a formula that ap- 
pears as Equation 7 in [11], which corrects a misprint in an earlier formula 
by Comtet and Fiolet [2]. 

The formula, as stated in Equation 25, uses some notation for (q + 1)- 
dimensional partitions. If (n, m) € N 9 x N is a nonzero tuple of nonnegative 
integers, then a (q + 1)- dimensional partition p of (n, m), denoted by p h 
(n, m), is a multiset {(si, ti), . . . , (s r , i r )} of nonzero tuples in N q x N that 
sum to (n, m) when counting multiplicities. We let \p\ = r denote the 
number of terms in the partition p, counting the multiplicity n p - s ,t of each 
tuple (s, t) in p. 
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Let y be implicitly defined by g as in Equation 1. We introduce the 
shorthands 

Q\n\ <9l s l+*£ 

Vn = y n (x) := ^r(x), 9s,t = <7s,t(x,n(x)) : = -0^yT^ y ^>- 

As the multiplicities fJ, p - s ,t sum to \p\ , the multinomial coefficient 
\p\ \ |p|! 



is well defined for any partition 

p = {(si,ti),...,(si,ti), ••• ,(si,^),---,(s«,t/)}. 

V v ' V v ' 

Corollary 7. 7/y and 5 are sufficiently smooth and related by Equation 1, 
then, for any xGU and nonzero n £ N q , 



n! ' - W V • • . »r,s,t, ■ ■ ■) At n I sm 5o,i (x, y(x)) 



ph(n,|p|-l) 
(0,1)^P 

Here the product is understood to be of fi p - s ,t copies for every distinct 
element (s, t) of the multiset p. The self-referring nature of the summation 
makes it not directly obvious that there is only a finite number of partitions 
p of this form for any n. Given such a partition p, only a < |n| of its 
elements (s,t) satisfy s / 0. Since p does not have (0,1) as an element, 
each of the final coordinates of the b remaining elements of p is at least two. 
Then 26 < \p\ — 1 = a + b — 1 implies that p contains at most \p\ = a + b < 
2a — 1 < 2|n| — 1 elements. This bound guarantees that any partition p 
should sum to a tuple smaller than (n, 2|n| — 1), implying that, for given 
n, there is but a finite number of multisets p of nonzero tuples in N q x N 
satisfying p h (n, \p\ — 1) and (0, 1) ^ p. 

On the other hand, the existence of partitions of this form can be seen by 
taking simple examples. For example, for q = 2 and n = (1,0), (2,0), (1, 1) 
one finds partitions 



n = (1,0): {(1,0,0)} h (1,0,0), 

n = (2,0):{(2,0,0)}h(2,0,0), 

{(1,0,1), (1,0,0)} h (2,0,1), 
{(0,0,2), (1,0,0), (1,0,0)} h (2,0,2), 
(26) n= (1,1): {(1,1,0)} h (1,1,0), 

{(1,0,0), (0,1,1)}, {(0,1,0), (l,0,l)}h (1,1,1), 
{(1,0,0), (0,1,0), (0,0,2)} h (1,1,2). 
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For these n, the Corollary claims 

_ 51,0,0 
2/1,0 - - 



5o,o,i 

52,o,o . n 5i,o,i5i,o,o 5o,o,25i,o,o 

52,o = H 2 2 3 ' 

5o,o,i 5o,o,i 35,0,1 

/ »n 5i,i,o . 5i,o,o5o,i,i . 5o, i,o5i,o,i 5i,o,o5o,i,o5o,o,2 

y^i) 5i,i= 1 2 1 2 3 • 

50,0,1 35,0,1 5o,o,i 1/0,0,1 

Coalescing the grid to a single point xo in Equation 23", one finds that 

5n(Xo) 

n! 



(28) ^ = E E 

0=po<pi<-<P|„|=n r'6T'(po,Pi,-.P|n| 




1 5sV*(xq,5(xq)) \ 
5 o,i(xo,y(x )) J ' 



where each star * is of type (s*,t*). Clearly Equations 25 and 28 are in 
a similar form. The only difference seems to be that in Equation 25 equal 
terms are grouped together into one term with a coefficient. It is not sur- 
prising that there are duplicate terms in Equation 28, as each term depends 
only on the types of the stars, not on how these stars are connected to form 
a tree. 

We first provide an example that introduces the flavor of the proof of 
Corollary 7. 

Example 8. Let n = (1, 1). Coalescing the grid to a single point (xq,Xq) 
in Equation 17, one finds 

/ ori N 1 5o,o,25i,o,o5o,i,o . 5i,o, i5o, i,o 91,1,0 

i^yj 51,1 = - o 3 1 o 

2 0o,o,i 90,0,1 9o,o,i 

_ 1 9o,o,25o,i,o5i,o,o 9o,i,i5i,o,o 
2 9o,o,i 9o,o,i 

at this point (xq,Xq). Clearly Equations 27 and 29 are equivalent. In this 
example we hint at how the terms in these equations are related, suggesting 
a link that generalizes to the generic construction in the proof of Corollary 
7. 

First of all note that in both Equations 27 and 29 the denominators can be 
determined from their numerators. Taking for granted that the coefficients 
agree, it therefore suffices to check that, for these equations, the monomials 
of the numerators of their terms agree. 

For each monomial in Equation 29, the orders of the derivatives in the 
numerators form a multiset p of triples in N 2 x N with (0,0,1) ^ p and 
p h (1, 1, \p\ — 1). It follows that every monomial in Equation 29 appears in 
Equation 27 as well. 

Conversely, we show that every monomial in Equation 27 appears in Equa- 
tion 29 as well, by pointing out which paths (0, 0) = po < pi < P2 = (1 5 1) 
and trees r' G T'(po, Pi, P2) correspond to it. For a given monomial in 
Equation 27, let p be the corresponding partition in Equation 26. To each 
triple (si,S2,t) in p we associate a star of type (si,S2,t). As can be seen 
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V 


{(0,0, 2), (1,0,0), (0,1,0)} 


{(0,1,1), (1,0,0)} 


{(1,0,1), (0,1,0)} 


{(1,1,0)} 


* 


(v+"i)(v + 

VI I 




(■,-+e 2 ) (-,.+6!) 

V I 


( v +ei) (v+e,) 

V 1 


(.,. + e 1 )(-,-+e 2 ) 

V 


t' 


r[ 

(0 !ei )( ei! n) 

o o 

V 

(O.ej.n) 


(0,e 2 )(e ; 

o 

y 

(0,62,1 


,n) 

o 

• 




'"a 

( e 2 > n ) 

o 

(0,e 2 ) 

V 

(0,e 2 ,n) 


( e l > n ) 

O 

(o, ei ) 

V 

(0 >ei ,n) 


(0,%) (ei,n) 

V 

(0, ei ,n) 



Table 1. For n = (1,1), the first row lists the multisets 
p with (0,0,1) ^ p and p h (n, |p| — 1). The second row 
depicts the stars ★ associated to each of these multisets. The 
third row shows the different trees t' that can be formed 
by connecting these stars, together with the labels of their 
vertices. 



in Table 1, one can, for each multiset p, connect these stars to a plane tree 
t' with two leaves, and sometimes there are several ways to do this. After 
demanding the first coordinate of the label of the left leaf to be (0, 0), there 
is only one way to label the leaves of r' that agrees with the types of stars 
of t' . Thus we find three trees t' 1: t' a ,t§ € T'((0, 0), (1, 0), (1, 1)) and two 
trees t' 2 ,t'^ 6 T' ((0,0), (0, 1), (1, 1)), each of which corresponds to a term in 
Equation 29. 

of Corollary 7. In both Equation 25 and 28, each term comprises some co- 
efficient and a monomial in the symbols g s ^t divided by a power of 50,1 of 
the same total degree. As this monomial uniquely determines the denom- 
inator, the following three steps suffice to show that Equations 25 and 28 
are equivalent. 

1. Every term in Equation 28 appears also in Equation 25. Consider an 
arbitrary term T in the right hand side of Equation 28. This term arises 
from picking a path = po < • • • < P| n | = n an d a tree t' G 7~'(po, • • • , P| n |)- 
Let p = {(si,ii), . . . , (sip|,i|p|)} be the multiset of types of the stars in r' . 
That is, p is the multiset of orders of the derivatives in the numerator of T. 

Since for any tree in 7~'(po, • • • , P|n|) t ne steps made by its leaves sum to 
n, it follows that s\-\ — • + S | P | = n - Moreover, with the exception of the root 
of t', the root of each star in Stars(r') connects to one of the t% + • • • + t| p | 
nonleaf descendants of the stars, implying that *! + ••■ + = \p\ — 1. As 
t 1 has no vertices with precisely one nonleaf descendant, none of the types 
in p can be equal to (0, 1). The orders of the derivatives in the numerator 
of Equation 28 therefore constitute a multiset p with elements in N 9 x N for 
which (0, 1) ^ p and p h (n, \p\ — 1). We conclude that, up to coefficients, 
each term in Equation 28 appears as a term in Equation 25 as well. 
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2. Every term in Equation 25 appears also in Equation 28. Suppose we 
are given a multiset p = {(si, t), . . . , (s\ p \, t\ p \)} of tuples in N 9 x N satisfying 
(0, 1) ^ p and p h (n, \p\ — 1) as in Equation 25. To each element (s, t) of p, we 
associate a star of type (s, t) whose labels are yet to be determined. Because 
the number of stars \p\ is one more than the sum of the nonleaf descendants 
ti + • • • + tu, one can always connect these stars to a plane tree r' with 
|si| + • • • + |s| p || = |n| leaves. Moreover, since each star in Stars(r') is of 
some type unequal to (0, 1), there is a unique path = po < • • • < P| n | = n 
and corresponding labeling of the leaves of r' such that r' € T'(po, • • • , P| n |)- 
We conclude that, up to coefficients, each term in Equation 25 appears as a 
term in Equation 28 as well. 

3. Corresponding terms have equal coefficients. Now we have shown that 
the monomials in Equation 25 are the same as those in Equation 28, it 
remains to show that their coefficients agree. Every term in Equation 28 
corresponding to a tree t' £ T'(po, ■ ■ ■ , Plnl); with p the multiset of the 
types of the stars of r', will contribute (— l)# Stars ( T ') = (— to the term 
in Equation 25 corresponding to p. The coefficients can be shown to agree, 
therefore, by counting, for every multiset p with (0, 1) ^ p h (n, \p\ — 1), the 
number of different plane trees that can be formed by connecting the stars 
of types corresponding to the elements of p. 

Let p be as in Equation 25. Let us call two elements (sj,tj), {sj,tj) of p 
equivalent whenever t i = tj . The equivalence classes form a new multiset p' 
in which each element [(sj,i)], or simply t for short, has multiplicity [it := 
Ss>o^p;s,<- Clearly p' has the same number of elements as p. Associate 
with each t £ p' a star with t descendants. Here we think of a leaf as a star 
with descendants. By Lemma 6, one can construct 

if b'l ) 

b'l Vw),mi, • • •/ 

different plane trees from these stars. For the vertices with t descendants of 
any such plane tree, the multinomial coefficient of {^ p - s ,t '■ s > 0} gives the 
number of ways to reinsert the leaves. It follows that we can form 

-( M )u( ) = L ( lpl ) 

\p'\ Vmo,Mi, • • ■) V- • • >Mp;s,t, • • ■/ bl V- • ■ ,V P ;s,t, ■ ■ ■) 

different plane trees from the stars corresponding to p. This agrees with 
the coefficient in Equation 25. We conclude that Equations 25 and 28 are 
equivalent. □ 
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